Determination of nucleation rates near the critical point 
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Abstract 

The nucleation rates derived for the condensation from a supersaturated 
vapor are examined both in the classical theory and in the modern coarse- 
grained field theory. By virtue of the scaling variable Xz it is shown that 
the method of steepest descent is irrelevant to evaluate the nucleation rate 
in the proximity of the critical point in the capillary approximation. If the 
logarithmic corrections to the activation energy of a droplet are taken into 
account, then the calculated nucleation rates provide an adequate description 

of the liquid-gas phase transition both near and out of the critical range. 
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The problem of homogeneous nucleation has been intensively investigated both theo- 
retically and experimentally (see, e.g., and references therein). The first theoretical 
approach to it, which is often referred to as classical nucleation theory, was worked out 
about 50 years ago mainly by Becker and Doring and Zeldovich ||. Later on the classical 
nucleation theory was generalized to a system of arbitrarily many degrees of freedom by 
Landauer and Swanson || and to the field theories by Calm and Hilliard |7j and Langer 
fill]. Langer pointed out that the coarse-grained procedure can be efficiently performed near 
the critical point where the radius of critical fluctuations is larger than the characteristic 
correlation length. 

Both classical and modern coarse-grained field theory assumes that the decay of an 
initially homogeneous metastable state should proceed via the formation of nucleated clusters 
of a new, stable state. The rate of relaxation of the metastable state, or nucleation rate, is 
given by the formula 

I = I exp (-/3AF e ) . (1) 

Here Iq is the preexponential factor, f3 = (ksT)' 1 , ks is Boltzmann's constant, T is the 
temperature of the system, and AF C is the excess free energy of the critically large cluster 
in the system. Even though the prefactors in both theories have different expressions, the 



numerical results for condensation from the supersaturated vapor [10 1 are very similar. In 



principle, it is possible to show [fnj that under the certain assumptions the prefactor derived 
in the classical theory may be obtained in a form equivalent to that of the field theory. To 
test the predictions of the nucleation theories a number of experiments on the separation 



of binary fluids near the critical point [12,13| have been performed. It turned out that the 



critical systems were more stable than it was expected from the theoretical calculations. To 



explain this difference between theory and experiment Binder and Stauffer |]I1,I5| argued 
that since the experimentalists measured the completion time of the transition, one has to 
consider both droplet formation and droplet evolution, growth and shrinkage, during the 
relaxation process. Near the critical point droplets grow very slowly and this circumstance 



decelerates the completion of the phase transition. A detailed description of the nucleation 
kinetics in near-critical fluids similar to the theory of coagulation by Lifshitz and Slyozov 
T6| , [iT| has been developed by Langer and Schwartz |T8| . The overall reaction rate was found 



to be much lower than the nucleation rate alone. 

In the present paper we would like to show that the nucleation rate itself calculated in the 
capillary approximation near the critical point is overestimated. Therefore, the agreement 
between theory and experiment should be even better. For the sake of simplicity, we will 
consider first the application of the classical nucleation theory to a first-order phase transition 
with only one order parameter, for instance, condensation from the supersaturated vapor. 

The excess of the Helmholtz free energy due to the formation of a spherical droplet of 



radius R is equal in the capillary (or thin wall) approximation [17j to the sum of the bulk 
and the surface energies 

Arc 

AF(R) = R 3 Ap + AnR 2 a , (2) 

3 

containing the difference in pressures inside and outside the droplet Ap and the surface 

tension a. Minimization of AF with respect to the radius yields the energy of the critically 

4 

large droplet AF C = -nR^a, where the critical radius is given by Laplace's formula R c = 



2a/ Ap. Using the reduced variables |T9[ Xz = R c y 47RX/3 and r = R/R c , we have 

(3AF = -\\\r* + A|r 2 (3) 

and 

(3AF C = |a| . (4) 

As first shown by Zeldovich ||, the size distribution function of droplets f(R) obeys the 
Fokker- Planck kinetic equation 

dt~ dR' [b) 

where 



J=-Bg + Af (6) 

is a so-called current in the size space, containing the diffusion and the drift coefficients B 
and A. The solution of Eq. ([5]) corresponding to J = Jq = is a well-known canonical 
distribution fo(R) oc exp (— (3AF). In terms of the ratio f(R)/fo(R), Eq. (|D may be 
rewritten as 



(7) 



The theory of first-order phase transitions is based on the requirement of the steady-state 
flux J = J ss = const. In the region of small r the distribution function f(R) is very close to 
the equilibrium distribution function fo(R), whereas it diminishes quickly for r > 1. These 
boundary conditions represent the constant replenishment of the stock of small droplets 
because of the thermodynamic fluctuations and the fact that droplets of supercritical sizes 
are removed from the system and merge into the new phase. The steady-state solution of 
Eq. ([7D that satisfies both boundary conditions reads 

J - 1= L bmr)- (8) 

In order to evaluate the integral analytically it is usually mentioned that the integrand has 
a very sharp maximum at R = R c due to the maximum of the activation droplet energy at 
this point. Therefore, one may replace the activation energy near the critical radius by its 
harmonic approximation 

1 / f) 2 /\F\ 

f3AF = p(AF) R=Rc + -p(-^r) (R-R c ) 2 (9) 

V / R=R C 

= f3AF c -\ 2 z {r-l) 2 

and then apply the method of steepest descent. We have 
B(r c )f (r c ) 



Js 



f3 fd 2 AF\ B(l)/ (1) Xz 



R c \| 2tt V Or 2 J r=rc R c v/? • 

Mathematically, the assumption that AF(r) has a sharp maximum at r c = 1 means that 



(10) 



P\AF»(1)\ 



■i -i 



« 1 , (11) 



where AF" denotes the second derivative of AF with respect to the reduced radius r. 
It is easy to see that in the capillary approximation 

/3|AF"(1)| = 2A|. (12) 

Our next step is to determine the dependence of the similarity number Xz on the critical 
exponents provided the temperature of the system approaches the critical one. 

As functions of 9 = 1 — T/T C , the thermodynamic quantities needed for our analysis have 
the following power-law approximations in the vicinity of the critical point: 

<7~0"' , Ap-O 13 ' 6 ' . (13) 

In classical theory the critical exponents (3' and 8' are equal to 1/2 and 3, respectively. The 
measured value of the exponent y! lies in the range 1.22—1.29 |2(| and we will use // = 1.25 
in our further calculations. Thus the critical radius tends to infinity at a critical temperature 
as = The similarity number near the critical temperature obeys the power 

law 

fSfj,' 



\ z = RjliraP ~ 9^ = 9^ (14) 



It is clear that the parameter Xz tends to zero at T — > T c and that the criterion (|il|) is not 
fulfilled. Therefore, in the thin wall approximation the method of steepest descent is not 
applicable to the evaluation of the nucleation rate near the critical point. 



In the semiphenomenological droplet model worked out by Fisher |2l| the activation free 
energy includes also the curvature term related to the small fluctuations in the shape of the 
droplet which do not change both the volume and the surface area of the droplet 

AF F (R) = --7rR 3 Ap + ATioR 2 + 3T/T 1 In — . (15) 
3 r 

Here r is the Fisher critical exponent, which is about 2.2, and ro is the radius of the smallest 

droplet in the system. It is worth noting that the spherical harmonic excitations of a droplet 



(Goldstone modes) have been calculated HQ also in the field theory. These calculations 



lead to the appearance of the logarithmic term in the expansion of the free-energy density, 
similar to the curvature term in Eq. ([1|), with the critical exponent r = 7/3. Then in the 
harmonic approximation 

/3AF F = /3AF[ - 9t + 2X z ( r _ 1)2 ; (16) 



/3AF/" = -r+^ + 3rlnS . (17) 
3 r 



Note that the critical radius R F appearing in Eq. (|T7|) is not the same as that given by 
Laplace's formula, but should be determined by the solution of cubic (with respect to R) 
equation. Now 

(3\[AF F (1)}"\ =9r + 2A| (18) 

and the method of steepest descent is relevant since the criterion flllD is fulfilled: 

(9r + 2A|)^ o -^(9r)- 1 «l. (19) 

Performing the saddle-point integration in Eq. (^|), one can find the nucleation rate in the 
droplet model approach 



T F_ B F (l)f F (l) 9r + 2A| 

Jm ~ — w c — V^^- (20) 

To compare the analytical expressions given by Eqs. ([10]) and ( p0|) with the numerical 
solutions of Eq. @, we plot in Fig. | the ratio ( J ss ) numeric / (</ ss ) ana i yt ic versus X z = R c \j4ira(3 
calculated both in the capillary approximation and in the Fisher droplet model. We see that 
formula ( pD| ) gives us values of the nucleation rate J ss that agree with the results of numerical 
calculations by Eq. ([8p within the 5% accuracy limit even for very small values of Xz- In 
contrast, in the capillary approximation presented by Eq. fllPf ), significant deviations from 
the numerical results start already at Xz = 2, which is assigned to a system rather far from 
the critical point. Thus the value Xz = 2 may be considered as a limit of applicability of the 

6 



classical expression (|T0D for the nucleation rate. At larger values of the similarity number 
both analytical expressions flIQ|) and (|20| ) fit well to the results of numerical calculations. 

The errors introduced by the method of steepest descent have been calculated numerically 
|23| for the condensation of a gas. It was found that the errors are negligibly small except 
for the smallest critical droplets and for the smallest values of the activation energy of the 
critical droplet. We show that this error can be parametrized by the single scaling variable 
Xz- Our analysis is valid for the Langer theory || also, in which the prefactor In is shown 
to be a product of the dynamical and statistical prefactors k and flo, respectively, 

J = {^fi . (21) 

The dynamical prefactor describes the exponential growth rate of the unstable mode at the 
saddle point. It is related to the single negative eigenvalue Ai of the generalized mobility 
matrix My = d 2 AF / d&d^j , where {^} is a set of macroscopic variables describing the 
system. The integration over a plane containing the saddle point is performed by the method 
of steepest descent and the statistical prefactor becomes 



fin = V 



2« \ 1/2 rdet(/3M /27r) 



,/3|Ai 



1/2 



(22) 



det(/3M'/27r) 

where V is the available phase-space volume of the saddle point, the index denotes the 
metastable state, and a prime indicates that the negative eigenvalue Ai as well as the zero 
eigenvalues of the matrix My is omitted. For the process of vapor condensation the simplified 
model contains only the order parameter. From the definition of the mobility matrix it 
follows that — (3Xi = 2A| [cf. Eq. fll2|)l. Therefore, Eq. ( p2|) cannot be applied in the capillary 



approximation for the liquid-gas system near the critical point because the criterion (Till ) is 
violated. To calculate the nucleation rate in the critical region one has either to evaluate the 



integral numerically or to insert the logarithmic corrections PJ22j| mentioned above in the 
activation energy of a droplet before the steepest-descent evaluation of the integral. In the 
latter case this curvature term will play a crucial role in the determination of the saddle point 
for the free-energy functional. This is the last important point in our discussion. Equation 



(H) can be rewritten also || via the imaginary part Im T of the analytic continuation of the 
free-energy density to the metastable state 

r 



7T 



hnF , (23) 



and the evaluation of ImjF has been widely discussed in the literature. Sometimes ImJ- 



has a form fl8| like x r exp(— x 2 ), where the dimensionless variable x 2 corresponds to the 
activation energy of a critical droplet times /3, and the saddle point is calculated in the 
capillary approximation. In the context of a steepest descent this implies that the term with 
the logarithmic corrections x T is considered as a slowly varying part of the integrand in the 
vicinity of saddle point. Therefore, these corrections are added to the activation energy after 
the saddle-point evaluation of the integral and the nucleation rate is overestimated again 
near the critical point. 
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FIG. 1. Ratio ( Js S ) numeric / (</ss) ana i yti c versus the parameter \ z = R c (47T0-/3) 7 corresponding 
to the Fisher droplet model (solid line) and to the capillary approximation (dashed line). 
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